Improvement of the 3ω\omega thermal conductivity measurement technique at nanoscale

The reduction of the thermal conductivity in nanostructures opens up the possibility of exploiting for thermoelectric purposes also materials such as silicon, which are cheap, available and sustainable but with a high thermal conductivity in their bulk form. The development of thermoelectric devices based on these innovative materials requires reliable techniques for the measurement of thermal conductivity on a nanometric scale. The approximations introduced by conventional techniques for thermal conductivity measurements can lead to unreliable results when applied to nanostructures, because heaters and temperature sensors needed for the measurement cannot have a negligible size, and therefore perturb the result. In this paper we focus on the 3$\omega$ technique, applied to the thermal conductivity measurement of suspended silicon nanomembranes. To overcome the approximations introduced by conventional analytical models used for the interpretation of the 3$\omega$ data, we propose to use a numerical solution, performed by means of finite element modeling, of the thermal and electrical transport equations. An excellent fit of the experimental data will be presented, discussed, and compared with an analytical model

Cryptic transcripts from a ubiquitous plasmid origin of replication confound tests for cis-regulatory function.

A vast amount of research on the regulation of gene expression has relied on plasmid reporter assays. In this study, we show that plasmids widely used for this purpose constitutively produce substantial amounts of RNA from a TATA-containing cryptic promoter within the origin of replication. Readthrough of these RNAs into the intended transcriptional unit potently stimulated reporter activity when the inserted test sequence contained a 3' splice site (ss). We show that two human sequences, originally reported to be internal ribosome entry sites and later to instead be promoters, mimic both types of element in dicistronic reporter assays by causing these cryptic readthrough transcripts to splice in patterns that allow efficient translation of the downstream cistron. Introduction of test sequences containing 3' ss into monocistronic luciferase reporter vectors widely used in the study of transcriptional regulation also created the false appearance of promoter function via the same mechanism. Across a large number of variants of these plasmids, we found a very highly significant correlation between reporter activity and levels of such spliced readthrough transcripts. Computational estimation of the frequency of cryptic 3' ss in genomic sequences suggests that misattribution of cis-regulatory function may be a common occurrence

A finite element method for fully nonlinear elliptic problems

We present a continuous finite element method for some examples of fully nonlinear elliptic equation. A key tool is the discretisation proposed in Lakkis & Pryer (2011, SISC) allowing us to work directly on the strong form of a linear PDE. An added benefit to making use of this discretisation method is that a recovered (finite element) Hessian is a biproduct of the solution process. We build on the linear basis and ultimately construct two different methodologies for the solution of second order fully nonlinear PDEs. Benchmark numerical results illustrate the convergence properties of the scheme for some test problems including the Monge-Amp\`ere equation and Pucci's equation.Comment: 22 pages, 31 figure

Topological Optimization of the Evaluation of Finite Element Matrices

We present a topological framework for finding low-flop algorithms for evaluating element stiffness matrices associated with multilinear forms for finite element methods posed over straight-sided affine domains. This framework relies on phrasing the computation on each element as the contraction of each collection of reference element tensors with an element-specific geometric tensor. We then present a new concept of complexity-reducing relations that serve as distance relations between these reference element tensors. This notion sets up a graph-theoretic context in which we may find an optimized algorithm by computing a minimum spanning tree. We present experimental results for some common multilinear forms showing significant reductions in operation count and also discuss some efficient algorithms for building the graph we use for the optimization

Pressure-induced and Composition-induced Structural Quantum Phase Transition in the Cubic Superconductor (Sr/Ca)_3Ir_4Sn_{13}

We show that the quasi-skutterudite superconductor Sr_3Ir_4Sn_{13} undergoes a structural transition from a simple cubic parent structure, the I-phase, to a superlattice variant, the I'-phase, which has a lattice parameter twice that of the high temperature phase. We argue that the superlattice distortion is associated with a charge density wave transition of the conduction electron system and demonstrate that the superlattice transition temperature T* can be suppressed to zero by combining chemical and physical pressure. This enables the first comprehensive investigation of a superlattice quantum phase transition and its interplay with superconductivity in a cubic charge density wave system.Comment: 4 figures, 5 pages (excluding supplementary material). To be published in Phys. Rev. Let

Automated derivation of the adjoint of high-level transient finite element programs

In this paper we demonstrate a new technique for deriving discrete adjoint and tangent linear models of finite element models. The technique is significantly more efficient and automatic than standard algorithmic differentiation techniques. The approach relies on a high-level symbolic representation of the forward problem. In contrast to developing a model directly in Fortran or C++, high-level systems allow the developer to express the variational problems to be solved in near-mathematical notation. As such, these systems have a key advantage: since the mathematical structure of the problem is preserved, they are more amenable to automated analysis and manipulation. The framework introduced here is implemented in a freely available software package named dolfin-adjoint, based on the FEniCS Project. Our approach to automated adjoint derivation relies on run-time annotation of the temporal structure of the model, and employs the FEniCS finite element form compiler to automatically generate the low-level code for the derived models. The approach requires only trivial changes to a large class of forward models, including complicated time-dependent nonlinear models. The adjoint model automatically employs optimal checkpointing schemes to mitigate storage requirements for nonlinear models, without any user management or intervention. Furthermore, both the tangent linear and adjoint models naturally work in parallel, without any need to differentiate through calls to MPI or to parse OpenMP directives. The generality, applicability and efficiency of the approach are demonstrated with examples from a wide range of scientific applications

Dynamics of skyrmionic states in confined helimagnetic nanostructures

In confined helimagnetic nanostructures, skyrmionic states in the form of incomplete and isolated skyrmion states can emerge as the ground state in absence of both external magnetic field and magnetocrystalline anisotropy. In this work, we study the dynamic properties (resonance frequencies and corresponding eigenmodes) of skyrmionic states in thin film FeGe disk samples. We employ two different methods in finite-element based micromagnetic simulation: eigenvalue and ringdown method. The eigenvalue method allows us to identify all resonance frequencies and corresponding eigenmodes that can exist in the simulated system. However, using a particular experimentally feasible excitation can excite only a limited set of eigenmodes. Because of that, we perform ringdown simulations that resemble the experimental setup using both in-plane and out-of-plane excitations. In addition, we report the nonlinear dependence of resonance frequencies on the external magnetic bias field and disk sample diameter and discuss the possible reversal mode of skyrmionic states. We compare the power spectral densities of incomplete skyrmion and isolated skyrmion states and observe several key differences that can contribute to the experimental identification of the state present in the sample. We measure the FeGe Gilbert damping, and using its value we determine what eigenmodes can be expected to be observed in experiments. Finally, we show that neglecting the demagnetization energy contribution or ignoring the magnetization variation in the out-of-film direction—although not changing the eigenmode's magnetization dynamics significantly—changes their resonance frequencies substantially. Apart from contributing to the understanding of skyrmionic states physics, this systematic work can be used as a guide for the experimental identification of skyrmionic states in confined helimagnetic nanostructures

A computational model of open-irrigated radiofrequency catheter ablation accounting for mechanical properties of the cardiac tissue

Radiofrequency catheter ablation (RFCA) is an effective treatment for cardiac arrhythmias. Although generally safe, it is not completely exempt from the risk of complications. The great flexibility of computational models can be a major asset in optimizing interventional strategies, if they can produce sufficiently precise estimations of the generated lesion for a given ablation protocol. This requires an accurate description of the catheter tip and the cardiac tissue. In particular, the deformation of the tissue under the catheter pressure during the ablation is an important aspect that is overlooked in the existing literature, that resorts to a sharp insertion of the catheter into an undeformed geometry. As the lesion size depends on the power dissipated in the tissue, and the latter depends on the percentage of the electrode surface in contact with the tissue itself, the sharp insertion geometry has the tendency to overestimate the lesion obtained, especially when a larger force is applied to the catheter. In this paper we introduce a full 3D computational model that takes into account the tissue elasticity, and is able to capture the tissue deformation and realistic power dissipation in the tissue. Numerical results in FEniCS-HPC are provided to validate the model against experimental data, and to compare the lesions obtained with the new model and with the classical ones featuring a sharp electrode insertion in the tissue.La Caixa 2016 PhD grant to M. Leoni, and Abbott non-conditional grant to J.M. Guerra Ramo